SPhT/92-042 



The 0{n) model on a random surface: critical 
points and large order behaviour 



B. Eynard and J. Zinn- Justin^ 
Service de Physique Theorique^ de Saclay 
F-91191 Gif-sur-Yvette Cedex, FRANCE 

In this article we report a prehminary investigation of the large limit of a generalized 
one-matrix model which represents an 0{n) symmetric model on a random lattice. The 
model on a regular lattice is known to be critical only for —2 < n < 2. This is the 
situation we shall discuss also here, using steepest descent. We first determine the critical 
and multicritical points, recovering in particular results previously obtained by Kostov. We 
then calculate the scaling behaviour in the critical region when the cosmological constant 
is close to its critical value. Like for the multi-matrix models, all critical points can 
be classified in terms of two relatively prime integers p, q. In the parametrization p = 
(2m + l)q ± /, mj integers such that < I < q, the string susceptibility exponent is found 
to be 7string = —21/ {p + q — I)- When / = 1 we find that all results agree with those of the 
corresponding (p, q) string models, otherwise they are different. 

We finally explain how to derive the large order behaviour of the corresponding topo- 
logical expansion in the double scaling limit. 
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1. Introduction 

We discuss in ttiis article an 0{n) symmetric matrix model in the large limit. The 
model has a direct interpretation in terms of a gas of loops [0, characterized by an 0{n) 
index, drawn on Feynman diagrams. The partition function is initially given by an integral 
over n + 1 matrices, but the integral over n of them is gaussian. The gaussian integrations 
generate an unusual effective one-matrix model which, unlike the standard one-matrix 
model 0, cannot be solved by the standard orthogonal polynomial method 0. 

The resulting model is solved instead, in the spherical limit, by the method of steepest 
descent. The saddle point equations generalize similar equations for the usual one- matrix 
model, but the method of solution involves new technical considerations. We show in 
particular that when n is of the form n = —2 cos{np / q) , with p, q positive relatively prime 
integers, the trace of the resolvent of the matrix is the solution of an algebraic equation of 
degree q with polynomial coefficients. 

As expected from the analysis of two dimensional regular lattices, non-trivial critical 
behaviour is only found for —2 < n < 2. We exhibit a class of critical and multicritical 
points, recovering in particular results previously derived by Rostov We mainly discuss 
the case n = — 2 cos{7Tp/q), p, q being two relatively prime integers, although some consid- 
erations also apply to arbitrary values of n. When p has the form p = (2m + l)q±l we find 
results consistent with the {p,q) models (in the CFT classification [Q) as previously ob- 
tained in the multi-matrix models . This includes the unitary family. However in general 
the scaling properties are different (or at least the operator content). This is not totally 
surprising since in general these models cannot be directly related to multi-matrix models. 
There is however one exception, the two-matrix model representing the Ising model on a 
random triangular lattice. This model is identical to an 0(1) model with a special cubic 
potential. It will be investigated first in order to show that the results obtained by the 
method of orthogonal polynomials are recovered at leading order. 

We then discuss multicritical points in a general potential in the 0(1) model. We 
finally study the critical points of the general 0{n) models. 

For the one-matrix model the large order behaviour of the topological expansion has 
first been determined by a direct analysis of the differential equations satisfied by the 
scaling specific heat. It has been shown that perturbation series is divergent and for half of 
the models (this includes pure gravity) non-Borel summable. This property has later been 
related to the existence of non-trivial saddle points of the initial matrix integral in the 
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large limit (the equivalent of the instanton solution of quantum mechanics and field 
theory). In the non-Borel summable case these saddle points are real, a property related to 
the unbounded nature of the integrand in the scaling limit or equivalently to the instability 
of the corresponding statistical models by creation of surfaces of higher topologies. 

In the case of the multi-matrix model, the large order behaviour of the critical and 
tricritical Ising models has been completely determined, again by a direct analysis of the 
differential equations. For the general [p, q) model only the 2kl divergence of perturbation 
theory has been established. The application of the steepest descent method to the multi- 
matrix integrals instead is not easy, because the leading contributions cancel since now the 
square of the Vandermonde determinant is replaced by the product of two independent 
determinants involving the eigenvalues of the first and last matrices of the chain. 

We show that the models considered here, because they can be solved by steepest 
descent, allow for a direct determination of the large order behaviour. 

Since the discussion of the critical points involves many technical considerations, in 
this article we explain the methods, establish the scaling of the free energy, explicitly 
calculate the trace of the resolvent of the matrix in the scaling limit, and compare it to 
the resolvent of the (p, q) string models. In a next article we shall examine the critical 
behaviour in more details. 

2. The critical Ising model 

We first consider a special example of the two-matrix model with the partition function 
given by 

Z^/dM.dM.e-''"/*"^^'*^''^^^', (2.1) 

Ml, M2 being two N x N hermitian matrices and the potential U having the form 

U{Mi,M2) = ^{M^ + Ml)-cMiM2 + {M^ + Ml)/3, < c < 1. (2.2) 

This model can of course be exactly solved by the method of orthogonal polynomials 
We expect that it exhibits a double scaling limit representing the critical Ising model on a 
random surface. We solve it however here by a different method and only at leading order 
for N large. A comparison with the large limit of the exact result provides a check for 
the new method of solution which we can then apply to the general 0{n) model. 
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We first transform integral (pTTI), setting A = (Mi - M2)/2 and 5" = Mi + M2 + 1 + c. 
Substituting into tlie potential (|2.2|) we find 



trC7(Mi,M2) =tr [A^S + V{S)], 

with 

AV = S^/3 - 2cS'^ + (3c - 1)(1 + c)S + const.. (2.3) 

Eventually tlie matrix A will be replaced by n copies, the Ising model being identified with 
the 0(1) case of the 0{n) model. 

The gaussian integral over A can be performed and Z is now given by an integral over 
only one matrix S: 

Z = J d5[det(5®l)]-i/VWf)*'^^(^). (2.4) 
The integral over unitary matrices can still be performed, and the result is 

Z = |A^(A)ndA.e-W^)^(^^)n(^^ + ^^)"'^'' (2-5) 

i=l i,j 

where we denote by A (A) the usual Vandermonde determinant of the eigenvalues of S. In 
this form the model can no longer be solved by the method of orthogonal polynomials, 
but conversely it is much easier, as we shall show, to solve it, in the spherical limit, by the 
method of steepest descent and then to also find non-trivial saddle points. 

2.1. The saddle point equation 

Varying one eigenvalue of S we obtain the saddle point equation 

j-^^ J J J 

with 

4V'{X) = (A - 2cf - (c - 1)2 = (A - 3c + 1)(A - c - 1). (2.7) 

In addition to the repulsive potential between eigenvalues, already present in the one- 
matrix model, we have now an attractive and unbounded potential at —Xj. A potential 
barrier must separate the eigenvalues from their refiected images and thus a neighbourhood 
of the origin must be free of eigenvalues, all eigenvalues being for g positive on the positive 
real axis. 
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It is convenient to rewrite equation (|2.6| ) in terms of the trace ujo{z) of the resolvent 
of the matrix S, 

, , 1 1 1 1 , , 

In the large N limit the distribution of eigenvalues p(A) = A^~^ J2i ~ ^i) becomes a 



continuous function and equation (|2.6|) can be rewritten: 



uJo{z + iO) + uJo{z - iO) = -uJo{-z) + V'{z)/g. (2.9) 



The distribution p(A) is then given by 



p{X) = -^{uj{X + iO)-uo{X-iO)). 



A particular polynomial solution of the equation ( pT9| ) is 

u:,[z) = ^^{2V\z)-V\-z)). (2.10) 

Setting 

UJo{z) = UJr{z) + —Uj{z), (2.11) 

we then look for a one-cut solution oj^z) of the homogeneous equation which, as a conse- 
quence of equation ( p77|) and the definition of cuq, for z large behaves like 

oj{z) = -A{2V'{z) - V'{-z)) + {I2g)/z + O (z'^) ~ -z"^ . (2.12) 

A quadratic relation. We have now to solve the homogeneous equation 

uj{z + iQ)+uj{z-iQ) = -uj{-z), (2.13) 

with the hypothesis that uj{z) is a real analytic function with a unique cut on the real 
positive axis {z = Q excluded). 

This equation has an equivalent form. Let us consider the even function 

r{z) = \ [u'^iz) + uj'^{-z) + uj{z)uj{-z)\ , (2.14) 

and calculate its discontinuity across the cut on z > 0, where uj{—z) is regular: 

r{z + iO) - r{z - iO) = ^ [uj{z + iO) - uj{z - iO)] 

x[uj{z + iO)+Lj{z-iO)+Lj{-z)]. (2.15) 



The r.h.s. vanishes and thus r{z) is an even regular function, which taking into account 
the large z behaviour of uj{z) is a polynomial of degree 4. Moreover equation (2.15) is 
equivalent to equation ( |2.13| ). 

Direct derivation. Note that equation ( |2.14|) can be directly derived from equation 
( p.6|) by multiplying ( p.6| ) by l/{z — Xi) and summing over i. From the two identities 



A^2 ^ ^ _ A, A, - A, ''^ ' N 



N"^ ^ \z - Xi z + Xj Xi + X 



and neglecting a term of order we find for uJo{z): 

uoliz) + ul{-z) + uo{z)uq{-z) = ^ {V'{z)uJo{z) + V'{-z)uJo{-z)) + ri{z), 

where ri is a constant when V'{z) is a polynomial of degree 2: 

■V'{X)-V'{z) V'{X)-V'{-z) 



gr^{z)= / dAp(A) 



z — X z + X 



2c- / dA/9(A)A. 



Shifting ujq{z) we obtain an explicit expression for r{z). We note in particular that a 
variation of g translates for uj only into a variation of ri. Finally the same equation can 
also obtained from the loop equations (the equations of motion) in the large limit. 

2.2. General one-cut solution 

A cubic algebraic equation. Equation ( p.l3|) implies that the branch points are square 
root branch points because turning twice around them we return to the initial function. 
However in the second sheet we find also the cut oi uj{—z). If we pass through this new cut 
we arrive in general onto a new sheet and so on. Here, due to the special coefficients, we 
find only three different sheets, which implies that uj{z) satisfies an algebraic equation of 
third degree. This equation can immediately be obtained by multiplying equation ( |2.14| ) 
by uj{z) — uj{—z). We then find: 

uj^{z) - ?,r{z)uj{z) = uj^{-z) - ?,r{z)uj{~z) = 2s{z), (2.16) 

where s{z) is an even function which is everywhere regular because the first expression is 
regular for z < and the second for 2 > 0. It is therefore a polynomial. 
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General solution. The explicit solution of equation ( p.l6|) can be written 

oo{z) = e-2^'^/3u;+(2) +e2*^/3w_(^), (2.17a) 

u±{z) = {s{z) ± ^/K{z)f'\ (2.176) 

where A is the discriminant of the equation 

/^{z)=s^{z)-r^{z). (2.18) 

Because we look for solutions u){z) which are not even in z, we verify that we can write 

^ = z3-3/2(w(z) - w(-z))(u;2^z) +u;2(-z) + \u{z)uj{~z)). 



This expression in particular shows that v A is an odd function. It follows that uj- {z) — 
ujj^{—z) and therefore 

u±{z) = ±^ [eT2z^/3w(2) - e±'^"/3cc;(-^)l . (2.19) 
zv3 L J 

Our normalizations are such that uj± ~ z"^ . 
The reciprocal formulae are also useful 

r{z) - uj+{z)uj-{z), s{z) = \ {uj\{z) + uj^_{z)) . (2.20) 

It follows: 



/\ = \{ul{z)-ul{z)). (2.21) 



The expansion ( ^.12| ) of u}{z) for large z implies that vA behaves like z^ . For the general 
one-cut solution \/A must then have the form 



= 12V?>icz {z^ - e2) [{z^ - a^) {z^ - b^)] , 

where we have chosen 6 > a > 0. The condition that A(0) = implies the relation 

s^{0)=r^{0). (2.22) 

The determination of has been chosen such that 

u±{z-iO)= e^^'^'/^u^ (z + iO), for z G [a, 6] . 

6 



The relation uj+{z) = U-{—z) then automaticaUy imphes that uj{z) has in the first sheet 
a cut only along [a, 6]. Indeed 



Critical points. A critical point is generated by the confluence of two different zeros 
of A. Two situations can arise: 

(i) a = e, this is the case of an ordinary critical point of the one-matrix model, the 
determinant coming from the integration over A playing no special role. It just modifles the 
form of the potential. From the point of view of the Ising model this is the low temperature 
phase where all Ising degrees of freedom are frozen. 

(ii) a = 0, this is a new critical point speciflc to the structure of integral (|2.5|), 
and the only case we shall consider from now on. The condition a = implies that 

< tr >, which characterizes the Ising spin fluctuations, diverges. Indeed N < tr > 
is proportional to J2i j + ^j) and diverges only when some eigenvalue of S vanishes 
This argument is conflrmed in the continuum limit by a determination of the scaling 
properties of A < tr >. 

Finally the critical Ising model in the continuum limit is reached when both conflu- 
ences occur simultaneously, a = e = 0. Note, however, that in this limit the eigenvalue 
distribution approaches a singularity of the integrand. Therefore we can no longer be cer- 
tain that the steepest descent method is valid beyond the leading order for A large. We 
have made a few checks but this point requires a more detailed investigation. 

2. 3. Simple critical point 

General critical solutions. Solutions corresponding to these critical points will be 
characterized by the property that they have a cut only between the origin and z = b > 0. 
The functions uj± corresponding to these solutions have the following general form: 



where A and B are even polynomials. The functions r, s, A have then indeed the expected 
analytic properties: 



for 2: > uj{z-iO) =u+{z + iO) +U-{z + iO), 
for 2; < uj{z - iO) = uj{z + iO). 



u± = 




(2.23) 




(2.24a) 
(2.246) 



s{z) = 2z- 
\/A = ibz 



-4 



(^2 _ ^2)2 ^ g ^^2^ _ ^2^^ ^2 (^2 _ ^2^) _ ^4^3 



Vz^ - [b''B\3A -B) + {z^ - fe^) A\3B - A)] , 



(2.24c) 



7 



provided A{z) and B{z) vanish at z = 0. The degrees of the polynomials A, B are directly 
connected to the large z behaviour of uj± and thus to the degree of the potential V{S). 

The cubic potential. In the case of the potential (|2.3|) uj± ~ z"^ for ^ ^ oo. These 
conditions imply A{z) = z^ and -B(^) oc z^, and therefore the general solution ( |2.23| ) 
reduces to 

1/3 

= z^ [^J\ - h'^/z^ T ib/z^ {yi^l^±iba/zy (2.25) 

The criticality condition a = defines a line g = gdc) in the g,c plane, and thus the 
two parameters 6, a are functions of c. Note that the double zeros ±e of A are given by 
= 6^ (a — 1)^/ (3(7 — 1). They vanish for a = 1 which thus corresponds to the Ising model 
critical point, a case we examine later. We first discuss the generic case a 7^ 1. 

To relate the parameters 6, a with the parameters c, g which characterize the special 
potential (|2.3| ) of the Ising model we can compare the expansions of uj+{z) for z large up 
to order 1/z. Using equations ( ^.12| ) and ( pj.l7a| ) we find at the critical point b{3a — 1) = 
12V3c, 62(3(7 - 5)/9 = (3c - 1)(1 + c), and 

12g^V3 = b'(^^-^y (2.26) 

The conditions < c and g > imply 1/3 < a < 19/9. The limit cr ^ —00, which cannot 
be reached in the Ising model, corresponds to a quadratic potential V{S). For c = 0, 
i.e. o" = 1/3 we know the critical value of g from the solution of the one- matrix model, 
gc = l/(12-\/3), value which agrees with the result ( |2.26| ). 

Critical behaviour. The critical behaviour is obtained from the small z behaviour of 
uj±. For cr 7^ 1, 

uj{z) oc z^/\ 

For a = 1 the behaviour is different and is discussed in the section 2.5. The amplitude in 
front of z'^/^, which is needed below to determine the normalization of uj{z) in the scaling 
limit, depends on the determinations. We therefore choose to let z approach the origin 
on the negative imaginary axis: z = — zA, A ^ 0+, a choice we shall keep throughout the 
article. Then 

uj+i-iX) ~ -2^/3(1 - a)64/3A2/3, (2.27a) 
uj_ i-iX) ~ -2-i/3(a + 1)62/3a4/3. (2.276) 
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2.4- The scaling limit 

We now calculate uj{z) near the critical point, when the cosmological constant is close 
to its critical value, i.e. for 

\x = l- g/gd < 1 , 

in the scaling limit. Scaling functions are solutions of equation ( ^.13| ) with a cut for 
z > a > 0. The scaling functions uj± have the general form 

uJsc,± = (Va2 - z2 ^ (c{z) v/a2 - ^2 ± izD{z)^ , (2.28) 

where D, C are two even polynomials whose form is fixed by comparing the large z be- 
haviour of Use with the small z expansion of the critical functions. 

Equations ( |2.27| ) yield the large z behaviour of iVsc for A +oo and thus the degrees 
of the polynomials C, D. Here C and D must be equal constants and thus: 

2 /3 

u;sc,± = -2-1/36^/3(1 - a) (Va2 - ± iz^ . (2.29) 

It remains to obtain the relation between a and x. One can for this purpose directly 
calculate the variation of the polynomials r, s at leading order in x. However we use here 
a different argument which can easily be generalized to more complicate cases. 

An auxiliary function. We consider the function Q,{z): 

d 

niz) = —{guoiz)). (2.30) 

We note that ( |2.11| ) implies 

^_du{z)_ 
12 dg ■ 

It thus satisfies the homogeneous equation (|2.13|) . It is then convenient to introduce the 



3T2'"/30(2) - e±2*"/30(-2) 



corresponding two functions 0±: 

"-^^^ 

and therefore 

n{z) = e-2^"/30+(^) + e2^"/30_(^). (2.31) 

Since Q is the derivative of a function which has a singularity of the form {z — ^o)^/^, 
zq = ±a, ±6, it can have a stronger singularity of {z — zq)~^I'^ type. Moreover from the 
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definition of uq we infer the behaviour of O for z large, fl{z) ~ 1/z. These conditions 
determine Q{z) uniquely as a function only of the location a, b of the singularities, a 
property we shall use later. At the critical point Q± have the form (|2.23| ), where however 
now A{z) can have poles at z = ±b and A{z) and B(z) no longer necessarily vanish at 
z = 0. The large z behaviour of Q{z) implies 

B{z) ~ 1/(673), A{z)=A{z)/{z^ -b"^) with A{z) = 0(l). 

Z — >CXD 

The leading correction to the large z behaviour of Wsc is found from equation (|2.29| ) to be 
of order z~'^^^. This determines the singularity at 2; = and implies A{0) = B{0). The 
unique solution is thus 



Therefore 



uJ±{z,g) = uj±{z,g,)TxgAiV3 (^1 - b^ / z^ ±ib/ z) + O [x^) . (2.33) 

z^l — b^ jz'^ ^ ' 

By identifying the large z expansion of Wsc with the small z expansion of (|2.33|) we find 



the relation between the parameter a and x. For example for z = —iX and A 0_|_ we 



derive from the expansion (|2.33|) : 

w+(-zA) ~ -2^/3(1 - (t)64/3a2/3(1 + 0{x)), 
u_{-iX) ^ -xg,?,^/^2^'%-^'^\-^/\ 

while the scaling function for A +00 

_ = -i(l - a)64/3a4/3A-2/3 + q {\-^'^) . 

It follows 

a^iV-2{^lI^^f\-'\ (2.35) 
Note finally that, taking into account the phase factors, ujsc can then also be written 

2/3 / , \ 2/3" 



(2.36) 



-2 + — a^j + \y—z — y z"^ — a? 
This expression can be conveniently parametrized setting 

2 = -acosv?, Wsc(2) = 22/^a^/^^>^/^(l-(7)cos(2(/?/3). (2.37) 
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We have assumed G [0, tt] for z G [—a, +a]. 

Interpretation: the resolvent in the (p, 3) model. We recall that pure gravity is a (2,3) 
model in the CFT classification, which can be constructed as the solution of the canonical 
commutation relation [P, Q] = 1 where P, Q are two differential operators of order 2, 3 
respectively [p!o| . 

In the corresponding one-matrix model the saddle point equations have a natural 
interpretation in terms of the resolvent of the operator of order two. However, as we show 
in section 4, the trace of the resolvent of the other differential operator is also solution 
of an algebraic equation, of third degree. Solving the equation explicitly, we obtain an 
expression proportional to expression ( p.36| ). Moreover we verify that the scaling relation 
( p.35|) is consistent with the behaviour expected for the specific heat u{x) in pure gravity. 

Our results are thus consistent with the hypothesis that the matrix S becomes, in the 
scaling limit, the differential operator of order 3, d"^ — (3/4) {w, d} of pure gravity. 

2.5. The critical Ising model 

The line of pure gravity reaches the Ising critical point when three zeros of A coincide: 
a = e = 0. This situation is found for a = 1 and thus c = (—1 + 2-\/7)/27, gc = 
5 X = lOc^ or 6 = GVSc. 

The function lo+ then becomes 



+ib/zY^\ (2.38) 



For z small for example z = —iX, X 0+, we find that co- gives the leading contribution 
to the two first terms of small x expansion: 



.(26)2/3^4/3 _ ^^^31/228/35-1/3^-2/3 ^ Q (^2) _ (2.39) 



This implies that u)±{z) in the scaling limit, g — gc small, have the form 

/ / \4/3 



Again ^0 s-nd a can be determined by expanding the scaling form for z large z = — iA, 
A +00: 



a;.c,-(.) = M2A)^/3 + ^) + ^ (^"'^' 
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Comparing with the expansion ( |2.39|) we conclude 

— _ /'?>/0^2/3 



a/h 



- my 
51/23-3/2^1/2 



(2.40a) 
(2.406) 



As we show in section 4 these results agree, up to the normalizations of x and w, with 
the scaling of the resolvent of the differential operator P = d"^ — (3/2){d, m(x)} of the 
(3, 4) model in the spherical limit. We thus find complete consistency with the established 
results of the critical Ising model in the spherical limit. 

Note that, here again, it is convenient to use the parametrization ( |2.37| ): 



(^—z + \/ z'^ — + ^— z — a/z^ — j 



4/3' 



and thus 



Z = -aCOSif, UJsciz) = -2^/^62/3(^4/3 (^4^/3^ 



^2.41) 



2. 6. The singular free energy 



We show now how one can explicitly calculate the singular part of the free energy 
F — InZ and thus verify its universal character, when compared to other quantities. 
We start from 



dF r 

9^-Q^= N{trViS)) ^N^ J dspis)Vis) 
f 

= — f dzuoiz)V{z) , 
2in J 



(2.42) 



where the contour in the last integral encloses the cut of ivoiz). 
We then multiply by g and differentiate again 



f{g) = ^fdzV{z)n{z) 



[2.43a) 
(2.436) 



where we have introduced the function Q{z) defined by equation ( |2.30| ). Using then the 
decomposition (|2.31| ) into ( |2.43| ) we find that we can write 



/ 



2mg 



dz 



-2i7r/3 



V{z)-e'^'^'^V{-z) \ n+{z). 



(2.44) 
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The singular part of the free energy is thus related to the singular part of Q{z) for a small 
which we determine now. 

Expansion of Q{z) for a? small. We recall that Vt{z) satisfies the homogeneous equa- 
tion ( p.l3| ). Moreover has a cut along [a, h] and behaves for large argument z like 1/ z. 
It has singularities of the form (a — z)~^^'^ and {z — 6)^/^ at z = ±a, ±6. It can be directly 
determined by solving an equation of the form (|2.16|) 

n^{z) -3r{z)n{z) ~2s{z) =0. (2.45) 

The functions r, s which appear here can have simple poles at 2; = a, b. This implies that 
they can be parametrized as 

, , 1 z^ + T^d , , 2 br^ 

^(^) = oT^ ^TT^ T^' ^(^) 



The discriminant A of the equation A = — r"^ has as a numerator a polynomial of degree 
six which must be a perfect square, condition which determines the parameters d and r. 
However, it is more transparent to use a different method. 
The critical function (a = 0) is given by equation ( |2.32| ): 

2/3 



The most singular, for z small, of the two functions behaves like z~'^^^. This determines 
the corresponding scaling function. Indeed it has the general form (|2.28| ) , C can have poles 
at z = ±a, and C, D behave like for z large. Thus 

Va — z^ ^ ^ 

with 

^0 = 2^/33-^26-1/3. 

Expanding Osc at next order for z large we get a correction of order z~'^^^. More precisely 
for A +00 

nsc,+(-a) ~ 3-'/'a^/'b-'/'X-^/'. 
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This next to leading contribution must also be the small z behaviour of the correction to 
the critical function we are looking for. This correction on the other hand is at most of 
order Xjz^ for z large. These conditions determine it entirely and we get 



z^Jl — jz^ ^ 



(2.46) 



Note that the expression in the r.h.s. has a singularity at z = 6 only of strength {^z — 
while {^z — 6)""^/^ could have been expected. Actually it can be verified that the location 
of the singularity at z = 6 varies by a term of order at least and thus the corresponding 
contribution to O is even smaller. 



The free energy. We note the identity 



dz 



3 ^vi-^V^ 



1 - 67^2 + ihjz 



-1/3 



Introducing the expansion (12^461) into equation (|]4|), using this identity and integrating 
by parts we find 



/(a) = /(O) ^ ^ ^ dz 



2z7r(7 



-1/3 



We now use the identity 



1 



consequence of the definitions (|2.10| ), (|2.11| ), in the integral. 

We need U- only at leading order. From ( |2.38| ) we deduce 



a;_ =z^ [^X-h'^lz'^ ^ihjz 



-2/3 



We then see that uj- gives no contribution to the integral. We calculate the contribution 
of a;o(^) by taking the residue at infinity. We need only the large z behaviour of io^iz) at 
leading order. 



ig\pi 



+ 0{z-^). 
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We thus obtain 

This completes the calculation of second derivative of the singular free energy Fs^ and we 
find: 

glF^^ig) ~ F"{x) ~ ^ j AzV{z)n,^ = -3 x 2-"^N\a/hf/\ (2.47) 

with X = 1 — g / Qc- In the case of the (3, 2) ordinary critical point this yields 

F"{x) = -3(19/81 - a/9)^/2(l - a)-^/^N^x^/^. 

In particular for c = we recover the expected result, i.e. twice the value of the corre- 
sponding one-matrix model. 

For the critical Ising model instead we find 



F"{x) = -{5/2y/^N'^x^/^. 



2.7. Large order behaviour 



We now look for a non minimal saddle point of the integrand ( p.5| ) . We know from our 
experience of the one-matrix model, that the second lowest saddle point is obtained when 
we move only one eigenvalue from the edge of the distribution p, to the nearest stationary 
position. The variation of the action E is then: 

where E is given by 

>^ 1 1 >^ 1 



dX \ g N ^ X~Xj N ^ X + Xj 



In the continuum limit 



dX \ g 

^ (2w(A) +w(-A)). (2.48) 



12^7 
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We are looking for a zero of d'E/dX in the interval (0, a) and thus Xf is small. We can thus 
replace uj{z) by its scaling form. We then find the variation 6F of the free energy in the 
scaling limit. For the generic critical point, using (|2.37| ) and (|2.35| ) we find 



ST, oc a^^^ oc x^/^, =^ SF oc exp ^— const, x 



5/4 



in agreement with the result found in pure gravity. For the Ising critical point, using ( |2.41J ) 
and ( |2.40£| ) we obtain 

ST, oc a"^/^ oc x^/*", ^ SF (X exp ^-const. x'^^^^ . 

in agreement with the known result for the Ising model. 

2. 8. Multicritical points 

With a general potential V{S) higher order critical points can be generated. If the 
polynomials A, B, in expression ( |2.23| ) behave like: A{z) oc B{z) oc ;22m+2 ^ small, 
io{z) oc ;2(6"i+2)/3 r^Yie minimal example corresponds to 

r{z) = s{z) = {z/bf^^^^. 

If in addition A{z) / B{z) 1 then uj{z) oc 2:(^"^+'*)/^. The minimal model is 

2 /3 

U± = (z/6)2-+2 (^y/l-byz^ ± ib/z^ , 

r{z) = s{z) = {z/bf"'+'^ - 2{z/bf'^+^. 

(Note the change of normalization compared to the m = case.) 

Scaling region. It will be convenient in what follows to classify the critical models 
when uj{z) behaves for z ^ like 

u{z) oc ^P/^ 

in terms of two relatively prime integers p, g, in analogy with the CFT classification (see 
section 4). For all the models considered up to now q = 3, while p belongs to one of the 
two sets p = 6m + 2, 6m + 4. 

To determine the functions in the scaling limit in terms of the deviation x from the 
critical cosmological constant we now use the same method as in the case of the simple 
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critical points m = 0. We first expand uj{z) at first order in x and then look for the leading 
term for z ^ 0. The two cases have to be discussed separately. 

First case p = 6m + 2. At leading order for z = —iX, A — > 0-|- we find 

u+{z) ~ (-l)'"+i2i/3(A/6)2'^+2/3. 

Therefore the scaling function cogc has the form (|2.28| ) with C, D polynomials of degree 
2m. 

It is convenient to normalize the potential in such a way that the variation 6lj at 
leading order in x and for z large is 

duj± ~ ^ixb/z. 



Then 6u± , which is independent of the potential and thus proportional to 0± of equation 

(P3), is 

= {Vr^ ± . (2.49) 



Expanding 6uj- for A — > 0-|- we find 

5(x;_(-zA)~-x22/362/3a-2/3. 

This behaviour of Su determines completely the scaling function. It can be expressed in 
terms of the function ^?(^): 



dt{t + iz)"'-^/^{t-iz)'^ . (2.50) 

-iz 

Calculating the integral we first verify that the function '^{z) has a form ( |2.28D : 

^(-2) = (^Va'^ - z'^ + iz^ ^ (C{z) A/a2 - ^2 + izD{z) 



For A —>■ +00 we find 

^i-iX) ~ (-l)'"S(m + 1, m + 2/3)(2A)2"^+2/3^ 

^2m+4/3(2A)-2/3 

m + 2/3 

with 

r(a)r(/j) 
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It follows 

u,c,±i^) = do^{±z), (2.51) 

with 

2-2m-l/3^-2m-2/3 

~ ~S(m + l,m + 2/3)' 

r , -|l/(2m+4/3) 

a/6 = 2 2^/3(m + 2/3)S(m + l,m + 2/3)a; 
Using equation (|2.47| ) we conclude that the singular part of the free energy then scales like 

F"(x)ocxi/(3-+2). 

As we show in section 4 all these results are consistent with the behaviour found in the 
(6m + 2, 3) string model. Moreover equation (2.48) yields then the expected scaling for 
the "instanton" action 6E oc a;i+i/(6m+4)_ 

Case p = 6m + 4. In the A ^ 0-(- limit u- gives the leading contributions: 

^_ = (-l)-+l22/3(A/6)2-+4/3 _ ^22/3(6/A)2/3 + Q (A'^/^' 

We introduce now the function i){z): 



dt{t + iz)'^+^/^{t-iz)"' . (2.52) 

-iz 

For A +00 we find: 

9-2/3 2m+2 

^(-iX) = (-l)-S(m + 1, m + 4/3)(2A)2-+4/3 + X-V^ + O ( A'^/^ ) . 

m + 1 V / 

In the scaling limit we conclude that uj{z) must have the form ( |2.51| ) with now 

UJ±{z) = 1?o^(T2:), 

with 

^-2m-4/32-2m-2/3 

~ ~S(m + l,m + 4/3)' 
^ = 2 [(m + l)S(m + 1, m + 4/3)a;]^/(^"'+^^ . 

It follows that the specific heat F"{x) oc x^/^^'^^^\ Again these results are consistent with 
a (6m + 4, 3) model (see section 4). Agreement is also found with the scaling form of the 
instanton action 6T, oc 
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3. The general 0{n) model 

We now generalize the previous model by replacing the matrix A by a set of n matrices 

d5dA...d^„e-''^/«'*^l^(-^? + - + -*"' + '^(^'l (3.1) 

V{S) being now also a general polynomial potential. 

The model has then an 0{n) symmetry. The quantity F = InZ can be interpreted 
as the free energy of a gas of loops, each indexed by an integer i, z = 1, ...,n, drawn on a 
random lattice of the form of a Feynman diagram. 

The corresponding model on regular lattices can become critical only for —2 < n < 2. 
We shall verify that here that the integral is even defined only for n < 2. It is thus 
convenient to set n = —2 cos 9. Although we could restrict ourselves to the interval < 
^ < TT, it is convenient for book-keeping purpose, to consider all positive values of 9. Note 
that the case n = 0,9 = {2m + l)n/2 reduces to the standard one-matrix model. For n = 1 
we can assign 9 = pn/ 3 to the multicritical models when lj{z) oc z^^^; the critical Ising 
model thus corresponds to 9 = 4tv/3. 

The integral over the matrices A^'s is still gaussian and can be performed. We can 
then parametrize S in terms of a unitary transformation and its eigenvalues A^. After the 
integration over unitary transformations, the integral ( |3.1| ) becomes: 

iN/9)ViX,) 



-iVE[A] 

dAe 



with the effective action E 

3.1. The saddle-point equation 

In the planar limit — > oo, Z can be calculated by the steepest descent method. The 
saddle point equation is: 
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We again introduce the density of eigenvalues p(A) = " ^0; Hilbert's 

transform, the trace of the resolvent: 

P(A) 



dA- 



'z- A' 



In the large limit p(A) becomes a continuous function and uq a function analytic except 
when z belongs to the spectrum of 5", that is with a cut on a segment [a, b] of the real 
positive axis. 



Equation (3^) may be written in term of ujq- 



uo{X + iO)+uJo{X-iO) + nuo{-X) = -V'{X), (AG [a, 6]). (3.3) 

9 

This linear equation has a polynomial solution: 

uAz) = -^j^{2V'iz)-nV'i-z)). (3.4) 

Note that the cases n = ±2 are special and must be examined separately. The function 
uj{z), defined by 

ujq = uJr + uj/g, (3.5) 
then satisfies the homogeneous equation: 

uj{X + iO)+uj{X-iO) + nuj{-X)=0. (3.6) 

Since Uoiz) behaves as 1/^ for ^ large, u!{z) has the large z expansion: 

uiz) = --^{2V\z) - nV\-z)) + ^- + (z'^) . (3.7) 
yz — n z 

A quadratic relation. Let us introduce the following function: 

r{z) = uj^{z) + uj^{-z) + nu{z)uj{-z). (3.8) 

We verify, as in the case of the Ising model, that the discontinuity on the cut of r{z) 
vanishes as a consequence of equation (|3.6| ): 

r{z + iO) - r{z - iO) = [uj{z + iO) - u{z - zO)] 

X {u{z + iO) + ijj{z - zO) + nu{-z)\ = . 

Therefore r is an even function, analytic in the whole complex plane. The behaviour of 
u) for z large, then implies that r is a polynomial. Again equation ( |3.8| ) can be directly 
derived from the saddle point equation or the loop equation. An expression of r(2;) in 
terms of the potential follows. 
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3.2. The general solution 

We introduce two aiixiliary functions: 



u;+{z) 



UJ-iZ] 



je/2 



2sin^ 

i 



uiz)-e'^/^u{-z) 



such that u!+{—z) = uj-{z). Then 



U!+{z)u!- (z) = r{z). 



Conversely u!{z) is given in terms of U- by 



u{z) = - e'^/^u+{z)+e-"'/''uj-{z) 



0/2, 



Equation (|3.6| ) then is equivalent to the simple relations 



u±{z - iO) = e^'^u^iz + iO), 



(3.9) 



(3.10) 



(3.11) 



which, as in Ising case, imply that uj{z) has cuts only on the positive axis. 

If we now choose 9 such that e**^^ = ±1, q being a positive integer, then u!{z) satisfies 
an algebraic equation of degree q. We have to examine the two signs separately. 

Case e*'^^ = 1. This implies 9 = Trp/q, where p is an even integer. The function s{z), 



s{z) = Uu^l + u'L) 



(3.12) 



has no discontinuity on the cut, and is analytic in the whole complex plane. It is therefore 
an even polynomial of a degree determined by the degree of the potential . We have thus 
the two following algebraic equations: 



(3.13) 



The solution of these equations is then = s ± v A, with 



(3.14) 



Note that is thus an odd function. The function uj is then given by equation (|3.10| ). 
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Case e*'^^ = — 1 . This implies 9 = np/q^ where p is now an odd integer. The expressions 
are quite similar but the role of s{z) and \fK are formally exchanged. It is now the function 

«W = 4K-^-). (3.15) 



which has no discontinuity on the cut, and is analytic in the whole complex plane. It is 
therefore an odd polynomial. We have thus the two following algebraic equations: 

ujj^uj-=r{z), — c<;i = 2zs(z) . (3.16) 

The solution of these equations is ujj^ = \/A ±is, with 



Here instead v A is even. The function uj is still given by equation (|3.10[ ) 



One-cut solution. We still have to determine the coefficients of the polynomials r, s, A. 
They can be found from the additional condition that uj has only one cut [a, b] on the 
positive real axis. We can see from ( 3.10|) that the singularities of u are the single roots 



of A. We demand that except for a and b (and —a,—b by parity), all the roots of A are 
double, in such a way that A can be written: 



iz^-a''){z^-b^)M^iz), 



where M{z) is an odd or even polynomial depending on the different cases. Due to the 
special form of the conditions ( |3.11[ ) all one-cut solutions uj± can be factorized: 

co± = n±{z) \±zA{z) + B{z)^{z'^ - a'^){z'^ - b^) 

where A and B are even functions, rational fractions in general because fl±{z) may have 
zeros, and the function O, which has only singularities at ±a, ±6, is a "minimal" solution 
of: 

n+{z) = n-{-z), n±{z-io) = -e^'^n^{z + io). 

This factorization property is a consequence of the algebraic equation satisfied by u!(z) and 
may not necessarily hold when 9/n is not rational. 
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3. 3. Critical points 

We again consider only critical points for which a = 0. The function uj{z) at a critical 
point has then a cut for < z < b. The general form of such a solution is analogous to 
the form ( p.23|) of the case n = 1: 



u;±{z) = (^^1-b^/z^ ±ib/z'^ ^^'^ A{z)^l-b^/z^±ibB{z)/z , (3.17) 

where I, q are relatively prime integers with < / < g and A, B are polynomials which can 
be chosen even without loss of generality. Indeed the situation A, B odd is equivalent to 
A, B even with the change I ^ q — I. 

One immediately verifies that r(z), s{z) and A(z) are polynomials of a form consistent 
with a one-cut solution provided A and B vanish at ^ = 0. 

A minimal realization of a critical point with polynomial potentials is: 



uj^[z) = Ti^z/bf-^^^ (Vl - 67^2 ± " ^ (3.18^) 

u±{z) = (z/6)2"^+2 (^^l-byz^±ib/zy~^^\ (3.186) 



Then in both cases 

a;+(z-zO) =e^^(^-'/^)w_(z + zO) 9 = n{l-l/q). (3.19) 

The two cases q — I even and odd correspond to the two situations e*^^ = ±1. 
For z ^ we find that 

u{z) oc zP/'^ 

where p = (2m + l)q — I in the case (a) and p = (2m + l)q + I in the case (6). Note 
that the values of p are such that m can also be defined as the integer part of p/2q since 
p/ {2q) — 1 < m < p/ {2q). Finally we see that for book-keeping purpose it is convenient to 
assign the angle 9 = Tvp/q to the critical point characterized by the integers {p, q). 

3.4- Scaling region 

We want now to derive uj{z) in the scaling region when the variable x = 1 ~ g/gc 
which characterizes the deviation of the cosmological constant from its critical value is 
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small. Functions which satisfy equation ( |3.6|) and are singular only at z = ±a and z = co 
have the general form: 



± iz 



-l/q 



C{z)^d^ -z'^±izD{z) 



(3.20) 



where again C and D are even polynomials. A comparison between the large z behaviour 
of Wsc and the small z behaviour of uj at the critical point yields the degrees of C and D. 
This determines them completely only for the minimal critical points m = 0. We then find 



± iz 



To obtain the relation between x and a and completely determine the form of the poly- 
nomials C, D for multicritical points (m > 0) we then proceed as in the g = 3 case and 
calculate the deviation from the critical form at leading order for x small. 

Deviation from the critical form at leading order. We now calculate the deviation 
from the critical form at leading order in the variable x = 1 — g/ Qc which characterizes the 
deviation of the cosmological constant from its critical value. We normalize the potential 
in such a way that for z large the variation 5uj± is: 

5uj± ~ ^ixhj z . 

As in the case of the Ising model we introduce the function 0(2): 

?L{z) = —{guJo{z)). 



Equation (|3.5|) then implies 



niz) 



duj{z) 
dg 



(3.21) 



It thus satisfies the homogeneous equation ( p.6|) . We also introduce the decomposition 



n±{z) = ± 



2sin^ 



n{z) 



n{-z) 



n{z) = - e^^/20+(2) +e-*^/2fi_(; 



(3.22a) 
(3.226) 



Prom the definition of ujq we infer the behaviour of Q for z large, fl{z) ~ Moreover, 
since O is the derivative of a function which has singularities of the form (z — zq)^^^, 
zq = ±a, ±6, it can have a stronger singularity of {z — zo)~^^^ type. These conditions 
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determine Q{z) uniquely as a function of the location of the singularities, as in the Ising 
case. The variation dco is proportional to Q at the critical point. At the critical point O 
has the form ( |3.17| ). To obtain its complete form we need its small z behaviour. This 
behaviour must be consistent with the leading correction to the large z behaviour of Usc- 
Since O is independent of the potential we can compare it to the form of Usc in the case 
(a) for m — 0. The leading correction to u>sc is then of order z'-/''~^. The unique solution 
is then 

With this information we can now explicitly calculate the scaling functions for all critical 
points. 

Scaling function. From the preceding analysis we conclude that for z large Wsc satisfies 

oosciz) - const. zP/" = O (xz^/'i-^^ . (3.23) 

It is easy to verify that this fixes the polynomials C, D. We now show that the scaling 
function Wgc can then be expressed in terms of the function 'd{z) given by the integral 
representation 

d{z) = dt (t + (t - iz)"^ , (3.24) 



the proof relying on a verification of condition (|3.23 ). Calculating the integral ( |3.24| ) we 



first verify that the function d{z) has a form consistent with expression (|3.20| ): 

d{z) = (\f(F^^ + iz^ ^'^ (C{z) \/a? - z'^ + izD{z)^ , 

where C, D are two even polynomials of degree 2m and ±l/q = p/q — 2m — 1. Moreover if 
we introduce the parametrization of case (a), p = (2m + l)q — /, we find 



^{I^Z)= I dt(t + Zz)(P-«-')/(2'?)(t-i;2)(P-9+0/(2g)_ 



It follows 

^{l, z) - -z) = (2z)P/V"('^-^)/(2'^)cTpg, 

where we have set 
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We then expand ^?(±^) for z large. 

Large z expansion. For z = —iX large we find 

^ ^ V ^ V ^ r(p/q + l) ^ ^ m + 1 

+ O (^XP/i-^^-^^ , (3.26a) 

(n\\2m-p/q 2p/q-2'm 

^(a) ~ ^ . (3.266) 

p/q — m 

Case (a) . If p = (2m + l)q — /, and thus p/q < 2m + 1, t?(— z) is asymptotically larger 
than the correction to i}{z). Moreover 2m — p/q = l/q — 1. The solution is then 

tOsc,±i^) = ^o^i±z). (3.27) 

Moreover comparing the expansion ( |3.26|) with the expansions of the critical functions uj± 
and 6u±: 

u+i-iX) ~ (-l)"^+i2^/«(A/6)^/^ 
6uj-{-iX) ~ -x2^-^/'^{X/by/''-\ 
we obtain the normalization constant t?o and the relation between a and a; 

Q^'^'-'^'' = 2(-+^-3^)/^[(p + ,-o/?]v, 



where we have used the definition (|3.25 



Case (b). If p = (2m + 1)^ + /, and thus p/q > 2m + 1, the correction to 'd{z) is 
asymptotically larger than the correction to 'i}{~z). Moreover p/q — 2m — 2 = l/q — l. The 
solution is then 

Usc,±iz) = ^o^iTz), (3.28) 

with 



pq 



pqX . 



Note that in the set of variables p, g, I the behaviour of a takes the same form in both 
cases. 
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3. 5. The singular free energy 

We can find the singular part of the free energy, using the same method as in the Ising 
model case. We have shown that 

where 0(z) is the function (|3.21|) . Using the decomposition ( 3.22fc|) we can rewrite equation 

(EH) 

The critical function for a = is 

l-l/q 



The scaling region. Let us now consider the case a ^ but small. For z small the 
function ( |3.31| ) behaves like z^l^~^ . This, together with the other properties, determines 
the scaling form of ^{z) 



Va2 - ^2 V ^ ' ; ' " sin(^/2) ■ 

Conversely, as in the g = 3 case, the next to leading term in the large z expansion of O 
provides the additional information needed to completely determine the first correction to 
the critical function for a small. This correction behaves like 2"^"^/"^, therefore 

0+(a) - 0+(a = 0) oc ^ Ul-h^z^ + ih/z^ 

z^^l — jz^ ^ 

The leading correction to Osc,+ (z) is: 
We thus have: 

0+(a) - O+(0) ~ (t) , ,^ , ( v/r^6W + 

sm(6'/2) V6/ z^^JX-lP' jz^ ^ 

The identity 
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allows to cast this expression into the form 



With this expression we can integrate by parts integral (|3.30|) . The second derivative of 
the singular part of the free energy is then given by 

X (^e'^/^V'{z)+e-''/^V'{-z)y 

We then use the identity 

e^^/V(z) +e-^^/V(-^) = 2z^7sin^ (^e'^^/^ ^o(^) - e^'/'c^o(-^)) 

+ 4sin^6'w_(2). 



We need U- only at leading order as given by expressions ( p.l8|) . We see that the con- 



tribution to the integral coming from uj- vanishes. The contribution due to ujq can be 
calculated by taking the residue at infinity. Then only the leading behaviour of ujq for z 
large is relevant: 

2z,sin^ {e-^'/^uo{z)~e^'/'uo{-z)) ^ ^'9sin9 cos{9/2) ^ 

In terms of the variable x = 1 — g / Qc we finally obtain: 

»W = = -^'(«/"(2 - n)2'-«/, Qy"\ 

If we set q = 3 and Z = 1 we find 3 x {a/lif'^^ , in agreement with the Ising model 

result. 

We derive from this expression the scaling of the free energy for all critical points: 

97 

F':^{x) (X ^ 7stri„g = -— (3.32) 

^ p + q — I 

3.6. Large order behaviour 

The calculation of section 2.7 can easily be generalized. The variation of the action 

is: 
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and we have: 

2 X ^ 1 n \ ^ 1 

~ Iv 2-^ X - X . AT 2^ X 4- X . 



= N (^^-^ - 2u;o(X) - ULOoi-X) 
V 9 

= -N [2a;(A) + m^(-A)] /g. 

In the critical region g close to gc, the argument A in the integral remains of order a. We 
conclude that 5T, scales like 

result consistent with the scaling of the free energy. Thus the variation 5F of the free 
energy has the form 

6F oc exp (^-const. , (3.33) 

result consistent with the 2k\ behaviour at large orders of the topological expansion also 
found in other matrix models. 



4. The resolvent in {p, q) string models 

We have solved the saddle point equations of the 0{n) model by transforming them 
into an algebraic equation for the trace of the resolvent of some matrix in the large N 
limit. It is therefore interesting, for comparison purpose, to discuss the form of the resol- 
vent in multimatrix models. The string equations of the (p, g) models can be generated 
by constructing a representation of the canonical commutation relations in terms of two 
differential operators P, Q of order p, q respectively (p, q are relatively prime and not both 
odd): 

P = d^ - (p/4){d^-^«(x)} + J]{^,(x),d^-2^}, 

Q = d« - {q/A){d'^-\u{x)} + Y.{w,{x)A''-^'}, 

where d means d/dx and {., .} means anticommutator. Moreover u{x) is the specific heat, 
the second derivative of the free energy. 

Let us consider the trace of the resolvent of one of these operators, for example Q: 

uj{z) = tr {z-Qy^. 
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In the spherical, i.e. the semiclassical hmit, because the non commutation between space 
and derivative can be neglected, uj{z) is given by 

uj(z x) = - f dx' f — (4 1) 

^ ' ' J J iiy)'^ - {q/2)u{x'){iy)'i-^ + 2Z,^^w,ix')izy)<i-^^ - z- ^ ' ' 

The integral over y just selects a residue. 

The scaling relations in the (p, q) model, in the same limit, are 

w^{x) oc u\x), u{x) ~ 

Thus uj{z, x) has the scaling form: 

ij{z,x) = z^l'^u . 

Moreover in the critical limit, i.e. for z large, it behaves like 

uj(z, x) oc zP/'? + xz^/«-^ + o f 



Note immediately that these scaling properties are consistent with those of the 0{n) model 
studied above if and only if / = 1, i.e. p has the special form p = (2m + 1)^ ± 1. 

We have shown that in the case of the 0{n) models oji^z^ satisfies, in the large limit, 
an algebraic equation of a special type. We show now that in this special class of (p, q) 
models the resolvent satisfies the same equations. 

It is well known that in the case of the one-matrix model, which corresponds to g = 2, 
io[_z^ satisfies a second degree equation. More precisely, setting p = 2m + 1, one finds 

, , m + 1 r^"^ dw' w'"^ 
u{z) = 



TT y/z + U' 

and thus 

uj''{z) = {z + u)Pl{z,u), 
where Pm is proportional to the polynomial part of the large z expansion of z'^{l-\-u/ z)~^/'^ 

Analogous properties for q > 2 are maybe less well-known. 

Operators of third order: q = 3. A short calculation shows that the trace uj{z) of the 
resolvent of the operator Q = z(d"^ — (3/4){d, w}) is, in the spherical and scaling limit, 
proportional, up to a rescaling of x, to the expression ( |2.51|) found in section2 2.4,2.5. The 
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result relies, in particular, on the property that p which is even and relatively prime with 
g = 3, is necessarily of the form p = 6m + 2 or p = 6m + 4. 

Generic situation. In the general case q > 4 the explicit functional form of the operator 
Q depends on the (p, q) models. In the semiclassical limit the coefficients Wi in equation 
( [4 .11) become proportional to with in general p-dependent coefficients. However, it is 
easy to verify, using for example the string actions ^ in the semiclassical limit, that when 
p = (2m + 1)^ ± 1, the semiclassical form of the operator Q is m-independent. If for 
convenience we change the normalizations u{x) —2u{x), z ^ 2i'^z, and then set in 
(0): 

y = 2'u}'!'^ cos 

we find for the denominator 

y'^ - qu{x)y'^~^ + 2z ^ 2 (u'^/^ cos(qv?) 

Taking the residue in y we obtain for uj{z) the expression 

f n}'!'^ sin Lp 
<^\z) oc / ax- 



z . 



ui/"^ sin((j(/?) ' 
where ip is solution of the equation 

u'^/'^ cos{qip) = z . (4.2) 

It follows 

1/9 / ^- \ 1/g 

z^ — ui 



2m^/^ sini/:' = (^z + \/z'^ - - 



2iu'^/'^ sm{qifi) = 2^/z'^ - ui. 
We now use the scaling relation x = ti(p+9-i)/2 a^fi get s = \J — u^. We find 



iV? _ _ c^i/g 



with a = w'^/^. In the two situations p = (2m + 1)^ ± 1 (thus a = a;i/(^"^+^) or a = 
^(j/(2m(j+2g-2)-j ^^iq integration yields an algebraic expression which, up to normalizations, 
agrees with the results ( pi.24D ,( P?77| ),( |3.28| ) obtained for the corresponding {p,q) critical 
points of the 0{n) model. Note finally that these special (p, q) models include the unitary 
family {q ±l,q). 
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5. Conclusion 

We have exhibited all critical points of the 0{n) model on a random lattice generated 
by the Feynman diagram technique, in the case n = ~2cos{np/q). We have calculated, in 
the spherical limit using steepest descent, the resolvent and the singular free energy in the 
scaling limit. In particular we have found that if we parametrize p as p = (2m + l)q ± /, 
< I < q, 7string = —"^l/ip + q — I). For / > 1 this result is surprising since it differs from 
what is found for the (p, q) critical points of multimatrix models: 7 = —2/{p + q — 1). One 
possible interpretation is that the operator which, in multimatrix models is coupled to the 
cosmological constant, is not present here*. 

We have finally characterized the large order behaviour of the topological expansion 
confirming the expected 2k\ behaviour already found in other matrix models. 

We have then shown that in the special case p = {2m + l)q ± 1 all the results we 
have obtained are identical to those found in the corresponding (p, q) string models by 
orthogonal polynomial techniques, and otherwise they are different. 

Moreover the techniques developed here allow to obtain a number of additional results 



which will be presented in a separate article [ 1 1 1 . In particular we shall study the effect of 



other relevant operators, a question which has some subtle aspects since, as the example of 
the Ising model reveals, negative powers of the matrix S have sometimes to be considered. 
We shall characterize more precisely the large order behaviour. Finally we shall exhibit 
solutions for generic values of n between —2 and +2, and discuss the n ±2 limit. 

Acknowledgement. While this work was in progress, Kostov communicated us a 
preprint [jl2| where the same model was discussed, though from a slightly different point 
of view. 



* We thank I. Kostov for this remark. 
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